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Abstract 



The boundary effects in the open Hubbard chain with boundary fields are studied. 
The boundary string solutions of the Bethe ansatz equations that give rise to a 
wave functions localized at the boundary and exponentially decreasing away from 
the boundary are provided. In particular, it is shown that the correct ground state 
of the model at half-filling contains the boundary strings. 
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In view of the experimental possibilities 0, interest in the theoretical study of 
transport properties for one-dimensional electron systems has been considerably en- 
hanced during recent years (see ref. and references therein). Thus, theoretical 
investigations of Tomonaga-Luttinger liquids (see |3| and references therein) in the 
presence of boundaries and potential scatterers are of crucial importance. A con- 
spicuous example of Tomonaga-Luttinger liquids with boundaries is provided by the 
one-dimensional Hubbard model with open boundaries and boundary fields (open 
Hubbard chain with boundary fields) @, 0, ^ 0. The Hamiltonian of this model is 
given by H = H^. a > + H hoimdaxy , where 

tf£X° = "EE (l&tfc+i, + tfWia) +ujr n n n jh (1) 

3=1 <r=U 3=1 
-^boundary = ~M n lT + n ll) - M n ^t + n L[)- (2) 

Here U is the coupling constant, hi is the boundary field at site I G {1, L}, ipj a (resp. 
ifjjo.) denotes the annihilation (resp. the creation) operator of an electron with spin 

° £ {1\ 1} a t she j G {1, 2, • • • , L}, and rij a = Vjo-Vv is the number operator. This 
model has been solved exactly by Bethe ansatz |4], |^, ||, [Tj, and its long-distance 
properties have been analyzed by using boundary conformal field theory |§, §, 0- 
Also, the boundary S matrices of the model have been obtained ||. 

In this letter, focusing on the complex solutions of the Bethe ansatz equations, 
we proceed to study the boundary effects in the open Hubbard chain with bound- 
ary fields. Originally, such kind of work has been done in the case of open XX Z 
chain with boundary magnetic fields || |I0|| . In refs. || [H|, the authors found 



boundary string solutions to the Bethe ansatz equations. The boundary strings rep- 
resent boundary bound states which correspond to the wave functions localized at 
the boundaries. In particular, it was shown that, in the gapfull regime with appro- 
priate conditions on the boundary magnetic fields, the ground state contains 



the boundary 1-string. 

We will show that, for the open Hubbard chain with boundary fields of type (0), 
there exist boundary bound states of the electrons corresponding to each boundary 
fields. If the boundary fields is sufficiently large (hi > 1 for I = 1, L) then the 
rapidity configuration in the ground state of this model contains boundary strings. 

Let us recall the Bethe ansatz solutions of the open Hubbard chain 0, 0, |], 0. 
The energy is represented as 

TV 

E = -2 cos kj, (3) 

3=1 

where the charge rapidities kj are the solutions of the Bethe ansatz equations 

M 

e^^^MMkjM) = ne(2( S infc j -A 5 )/c)e(2(sinfc j + A 5 )/c), (4) 

5=1 

M N 

J] e ((A 7 - A*)/c) e ((A 7 + A s )/c) = ]Je (2(A 7 - sin kj)/c) e (2(A 7 + sin %)/c) . (5) 

5=1 j=l 
5(^7) 
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Here N and M respectively denote the number of electrons and electrons with down 
spin, and 



c = U/2, (6) 
i _ h e ~ ix 

.(«) = £±|. (8, 

Note that, in this model, the solutions of the Bethe ansatz equations are restricted 
as Re(fcj), Re(A 7 ) > and kj,A y ^ 0. 

Putting k-j = — fcj,A_ 7 = — A 7 and k$ = 0,Aq = 0, we can rewrite the above 
Bethe ansatz equations into the more tractable forms, 

M 

e &WW)p(k j ,h 1 )P(k j ,h L ) = e(2wik j /c)- 1 JJ e (2(sinJfe i - A s )/c) , (9) 

S=-M 

M N 

e(A 7 /c)^ I] e((A 7 - A 5 )/c) = [] e (2(A 7 - sin k^/c) . (10) 

5=-Af j=-JV 
5(^7) 

Also, the energy is represented as 

AT 

£ = - ^ COS% + 1. (11) 
j=-N 

For the problems with periodic boundary conditions, we usually adopt the string 
hypothesis fll2| . The string hypothesis states that, in the thermodynamic limit, the 
set of solutions {kj, A 7 } of the Bethe ansatz equations can be split into three kinds 
of solutions, which are real k/s, combination of n A 7 's (A-strings of length n) and 
combination of n A 7 's and 2n k/s (fc-A-strings of length n) ||12|| . (We shall call these 
the bulk string solutions.) The A-strings and the A-strings can be interpreted as 
some kind of bound states jI2], [T3|] . 

For the problems with open boundary conditions, since the equations (|9|) and 
(|H]|) are very similar to the one for the model with periodic boundary conditions, it is 
reasonable to conjecture that the bulk strings make up the (sub)set of the solutions 
of the Bethe ansatz equations. 
Boundary string 

For the model with open boundary conditions, however, there may exist the 
complex solutions which are the different kind to the bulk strings. Note that, for 
the model with open boundary conditions, the energy (0) is the integrals of motion 
but the total momentum is not. Thus YljLi kj does not have to be real, in contrast 
to the case of the periodic boundary conditions. The rapidity configurations which 
contain pure imaginary roots (or roots whose real parts are n) and do not contain 
the complex conjugate of these pure imaginary roots are admissible. (We shall 
call this kind of solutions which correspond to the boundary fields ' boundary string 
solutions', m [L0|| ) 
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One particle systems 

As a warm-up exercise, we first consider the one electron system (N = 1). The 
wave function of the Bethe ansatz states is given by 

$ CT ( n ) = A a (k)e ikn - B a (k)e~ ikn (12) 

where a =| or j is the spin of an electron. The relation between A a {k) and B a (k) is 
obtained by investigating the boundary scattering at left or right boundaries. The 
rapidity k is determined by the consistency condition for the left and right scatterings 
(Bethe ansatz equation). The boundary scattering at the boundary n = 1 gives 



= (13) 



and at n = L gives 



B a (k) = \^ e J^ ik2(L+1) Mk). (14) 



The consistency condition for both left and right scatterings gives 



1 - h x e lk 1 - h L e tk ^(l+i) ^ 



1 — h\e~ tk 1 — h L e~ 



ik 



There are two type of solutions for (fL5f): one is the real solution which correspond 
to the free particle interacting at the boundary, the other is the pure imaginary 
solution which can be interpreted as the bound state at the boundary. When k 
is pure imaginary, the right hand side of decreases exponentially fast in the 



thermodynamic limit (we deal with only the case Im(fc) > 0). Therefore, if eq. (15 ) 
has pure imaginary solution, the rapidity k must take values of which the left hand 
side of ( |T5| ) vanishes. Let km = ixci) > kin — ^X(l) be the pure imaginary solutions 



such as 

l-e-xm/ix ~ e -2x ( i)(i+i) } 

1 -e~ x ^h L ~ e - 2x w {L+1] . (16) 

Note that the above equations have solution with X(i) > only when hi > 1 (/ = 
1,L). The solution k^ {k(i)) corresponds to the wavefunction localized at the site 
1: $ CT ~ e~ nx W (resp. site L: $ CT ~ e- {L ~ n)x ^ ) 
Ground state 

Now we investigate the ground state of the open Hubbard chain with boundary 
fields at the half filling^]. For the model with periodic boundary conditions, in the 
ground state, the rapidities {kj, A 7 } are all real and fill the Fermi seas. However, as 
examined above, we must consider the possibility of the existence of the boundary 
strings. Our task is, then, to compare the energy of the configuration with boundary 
strings to that of the configuration without the boundary strings. (In the remainder 

1 In this letter, the term half-filling means that the real charge rapidities kj fill the Fermi-sea 
between — tt and n. 
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of this letter, we treat only the case that the real roots {kj, A 7 } fill the Fermi-seas.) 
To calculate the effect of the existence of the boundary strings to the energy, we 
must investigate, in addition to the energy of the boundary strings, the shift of the 
densities of the real roots due to the existence of the boundary strings. 

We shall introduce the densities of roots and holes. The number of allowed 
solutions for the Bethe ansatz equations ([D and (H) in the intervals (k, k + dk) and 
(A, A + dA) are expressed as L[p(k) + p h (k)]dk and L[a(A) + a h (A)]dA. Here p(k) 
and a (A) are the densities of roots (filled solutions), and p h (k) and a h (A) are the 
densities of holes (unfilled solutions). 

If there exist two boundary strings X(i) — ln/ii, X{L) = ln/ix for each boundary 
fields (and there do not exist holes in the real roots), the densities p(k), cr(A) satisfy 
the following integral equations 

1 r B 

- + 2 cos k / dAa(A)K(2(sm k - A)) 

7T J-B 



P {k) 



1 

L 



2tt 



(£(k, hi) + £(k, h L )) - 2 cos kK(2 sin k) 



(17) 



o-(A) = 2 ( Q dkp(k)K(2A - 2 sin k) - [ dA'a(A')K(A - A' 

J-Q J-B 



1 

+ L 

1 



2c 



-K{- 



2cA 



-) + 



2c 



-K{- 



2cA 



c — 2sinh%(i) c — 2sinhx(i) c + 2sinhx(i) c + 2sinhx(i) 
2c 2cA , 2c „, 2cA 



2 sinh X(L) 



-K( 



+ 



2 sinh c + 2 sinh 



-K( 



^(A), 



l c + 2sinhx(L) 

(18) 



where £(h,h) = (h 2 - 1)/(1 - 2/icosA; + /i 2 ), iT(x) = c/[vr(x 2 + c 2 )]. Here Q and 
i? are charge and spin pseudo Fermi-momentum, respectively. We assume that, to 
examine the densities of order 0(1/L) in the half filled case, we can take B = tt and 
Q = oo. 

If we separate the densities into the part of order 0(L°) and 0(1/ L) as 



p(k) = Po (k) + ^ Pl (k) + O(^) 
a(A) = a0 (A) + y i (A) + O(^) 



(19) 



then po(k) and cr (A) are identical to the densities of fc's and A's, respectively, for 
the ground state with periodic boundary conditions. The ground state energy of 
the model with periodic boundary conditions is, then, given by 



periodic 



—L 



dkp (k) cos k. 



(20) 



The 0(1/ L) contribution pi(k) can be divided into three parts which reflect the 
geometry (open boundary conditions), the effect of the boundary fields, and the 
existence of the boundary strings; 



Pl(k) = Pgeom(fc) + Pfield(^) + Pstring(&)- 



(21) 
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We further separate p fic id(£:) to p ne id(&) = Pmd(k, h) + p &c i d (k, h L ) where p &e id(k, hi) 
is the contribution from the boundary field hi, and p s t T (k) to p str (k) = p s t T {k, hi) + 
Pstr{k, hi) where p s t T (k, hi) is the contribution from the boundary string correspond- 
ing to the boundary field hi. 
Here p gC om{k) is given by 

COS A; f°° £ -iw sin k-c\w\/ 2 
PgcomO) = dW 1 + e c\w\ 2 C ° S kK ( 2 Sil1 fc )> ( 22 ) 

and p fie id(&, h) is given by 

P(h h\ m^ik r°° p -2iw(sink-sink')-c\w\ 

Pfrf, h) = dk> dw- — . (23) 



2ty 2ix 2 J-tt J-oo cosh(cw 

The form of p s t T {k, h) is depend on the value of the boundary field h. 
The energy in the case that there are boundary strings, is given by 

dkpo(k) cos k — / dkp±(k) cosk — 2(coshxi + cosIixl)- (24) 

-7T J — 7T 

Then the energy difference between the configuration which contains boundary 
string and the configuration which does not contain the boundary string, corre- 
sponding to the boundary field hi, is given by 

AE(hi) = -2 cosh X i + e(hi), (25) 

where 

dkp s tr(k, hi) cos k. (26) 

If the strength of the boundary field hi is such that c < |2 sinh 1 , then 

K ( Vi cA ) K(2 sin k - 2A) 

V 2 sinh Xi-cJ v > 

2 sinh Xi ~ c 

and 



roc JA . r J\ Z Sill ft — ZiV 

p str (Mz) = 4ccosA;/ rfA ^i^p£^ ^ (27) 

J-oo 2 smh Yi — c 



,i n 1 Z* 71 ,, /*°° , a 4ccos 2 /c r _ 2cA . r ^ . . , . « 

e(hi) = -- dk dA K( )K(2smk - 2A). (28) 

v ; LJ-n J-oo 2sinhY, -c v 2sinhv, - c ; v 1 y 1 



oo 2 sinh X/ — c 2 sinh c y 

Noticing > (x is real), we see that e{hi) < 0. 

On the other hand, if |2 sinh xz | < c, then 



(i h\ cosk 

Pstr{k,hi) = - 

where 



c c 
G{- - sinhxi, fc) + G{- + sinhxi, fc) 



(29) 



f°° , e -»«"™ y -*M x /-oo i^(c(A — sin y) / x) 

G(x,y) = / dw — = — dA . ; . 30 

v yy 7-oo coshcw/2 2tt7-oo cosh(vrA/c) v 1 
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In this case, e(hi) is given by 



e(hi) = — — / dk 



! ,' " < ° / f — G (?- sinn Xi,^) + G'(^ + sinhx«,fc) 



(31) 



Again, noticing G(x,y) > (x > 0), we see that e{hi) < 0. 

We thus conclude that the configuration which contains boundary strings has 
lower energy than the configuration which does not contain the boundary strings. 
Therefore, if hi > 1 for I — 1, L, the ground state rapidity configuration of the open 
Hubbard chain with boundary fields at the half-filling contains the boundary strings 
for each boundary fields. 

The surface energy of the Hubbard chain with boundary fields 

-^surface = E S — -Aperiodic) (32) 

where E gr is the ground state energy of the Hamiltonian H and -^periodic is the ground 
state energy of the model with periodic boundary conditions, is given by 



^surface = ~ dk[p geom (k) + p fle ld(&) + Artr(&)] COS k 

J — oo 

-2(cosh xi + cosh xl)- (33) 

Conclusion 

Although, we did not yet performed to classify the all solutions of the Bethe 
ansatz equations (Q) and (|^), we have found that there exist the boundary string 
solutions in addition to the bulk strings. 

It was known that the open Hubbard chain with another type of the bound- 
ary fields -^boundary = — ^i( n iT — n u) ~ ^i( n LT — n Li) can be solved. The related 
Bethe ansatz equations for this Hamiltonian contain factors which depend both the 
coupling constant and boundary fields 0, 0]. The analysis of these Bethe ansatz 
equations will be appeared elsewhere. 

After this letter was written up, we find the preprint which relates to our results 



14j| . In their paper, boundary strings for the both charge and spin rapidities are 
studied. 
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